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In most Yang-Mills models the vacuum where magnetic monopoles condense coincides with that 
where center vortices percolate, thus it is not clear which of these two properties is most directly 
involved in producing confinement. It is pointed out that there is a class of 3D gauge models, 
which can be though of as duals of Q-state Potts models with Q < 1, where the magnetic monopole 
condensation is a necessary but not sufficient condition for percolation of center vortices. A set 
of numerical tests at Q = shows that there is a vacuum in which the magnetic monopole 
condensate does not yield confinement, in the sense that large Wilson loops obey a perimeter law. 
In such a vacuum the center vortices form a dilute gas of loops. At stronger coupling there is also 
a truly confining vacuum where both confining mechanisms are present. 
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1. Introduction 

Magnetic monopoles [[Ij] and center vortices [§] are widely believed to be the most important 
degrees of freedom for confinement in Yang Mills theories. 

Plausibility arguments suggest that magnetic monopole condensation implies a dual Meissner 
effect which pushes out of the vacuum the colour field and gives the well-known physical picture of 
confinement in terms of dual Abrikosov vortices describing the confining strings joining the quark 
sources. Considerable evidence for this dual Meissner effect has been accumulated on the lattice 
[^], including the definition of a disorder parameter demonstrating the condensation of magnetic 
monopole below the deconfining temperature [Q] 

Center vortices are string-like excitations formed out of the center of the gauge group which 
are expected to encode all the infrared physics of confinement [|5[ [6|]. When they percolate [f7|], 
produce a very efficient disordering mechanism which could lead to the area law decay of large 
Wilson loops. 

In most YM models the phase with magnetic monopole condensation coincides with that 
where center vortices percolate, thus it is not clear which of these two properties is most directly in- 
volved in producing confinement, or, to be more precise, the area law decay of large Wilson loops. 
There are many open, intertwined questions about the validity of these confinement mechanisms. 
In particular, is it possible to derive monopole condensation from percolation of center vortices or 
vice versa? Are both mechanisms necessary for confinement? are they also sufficient? 

In this talk I try to answer these questions by studying a particularly simple class of 3D gauge 
models, which can be though of as duals of Q-state Potts models. In these models the confining 
mechanisms can be easily identified in some specific geometric properties of the random graphs 
associated to the configurations of the Q-state Potts models. In particular it will be evident that the 
percolation of center vortices implies the condensation of magnetic monopoles. On the contrary, it 
is pointed out that when Q < 1 the magnetic monopole condensation is not necessarily associated 
to the percolation of center vortices: it is demonstrated through a numerical experiment that there 
is a vacuum state in which, although the magnetic monopoles condense, quark sources are not 
confined, because large Wilson loops decay exponentially with the perimeter instead of the area. 
It is also shown that in this theory there is a confining vacuum only when the magnetic monopole 
condensation is associated to the percolation of center vortices. 

The contents of this contribution are as follows. In the next Section the main properties of the 
Q-state Potts model and of its gauge dual are described with a particular emphasis on the definition 
of Wilson loop, which plays an essential role in the studies of confinement. In the following Section 
the nature of the vacua of these gauge duals in the range < Q < 1 and the salient features of the 
phase diagram are discussed. In Section 4 a local Monte Carlo algorithm for simulating Potts 
models in the range < Q < 1 is presented. Finally, in the last Section, some numerical results 
generated with this algorithm are reported. 

2. Q state Potts model and its dual 

The Hamiltonian of the (ferromagnetic) Q-state Potts model is, for Q integer, H = — L(r j) ^oj o> > 
where the site variable a, takes the values a, = 1 , 2, . . . , Q, with (ij) ranging over the links of an 
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arbitrary lattice or graph A. This model is symmetric under Sq, the group of permutations of Q el- 
ements. The canonical partition function Z = e~$ H can be rewritten in the Fortuin Kasteleyn 
(FK) random cluster representation: 

Z= £ w G = ^{b,c)v b Q c , (2.1) 

GCA b,c 

where v = e$ — 1 and the summation is over all spanning subgraphs of A, wq = v h Q c is their weight 
expressed in terms of the number b of edges of G, called bonds, and the number c of connected 
components, called FK clusters; £l(b,c) is their multiplicity. This representation now defines a 
model for any real or complex Q, which acts as the fugacity controlling the number of FK clusters. 

When A is a three-dimensional lattice, the FK random cluster representation is also useful to 
define a gauge dual of this spin model (for any complex Q) . The gauge dual lives in the dual 
lattice A. The most basic observables of any gauge model are the Wilson loops. In the present case 
these are associated to the closed paths 7 £ A. For any spanning subgraph G C A the Wilson loop 
Wy measures the topological entanglement between 7 and G. More precisely we attribute to the 
Wilson loop Wy(G) the value 1 if no FK cluster of G is linked to 7, otherwise we set Wy(G) = 0. 
The vacuum expectation value of Wy is defined accordingly: 

(Wy) = Y,Wy(G)Cl(b,c)v b Q c /Z . (2.2) 

b,c 

In the special cases where Q = 2,3, . . . this definition coincides with the one obtained by apply- 
ing the usual Kramers-Wannier duality, provided one defines the topological linking as a winding 
modulo Q [§]. The gauge theory dual to Q = 1 Potts model, corresponding to random percolation, 
has been studied in detail in In particular it has been shown that, although the gauge group 
is trivial, it behaves like a full-fledged gauge theory with a confining vacuum (corresponding to 
the percolating phase), a string tension having a well-behaved continuum limit, a non trivial glue- 
ball spectrum [ |i"0| ] and a deconfinement transition at a well determined temperature. In this talk I 
describe some new features of this kind of gauge models in the range < Q < 1. 

3. Confining vacua 

Wilson loops provide us with a fundamental tool for a precise definition of confinement in a 
pure gauge theory. A confining phase is expected to show up in an area law decay for the vacuum 
expectation value of large Wilson loops (Wy) . This exactly means that if 7 is scaled up keeping its 
shape fixed and increasing the area A of the encircled minimal surface, then (Wy) °< e~ aA , where 
a defines the string tension. 

For a generic Q > one expects two kinds of vacua, depending on the value of /3 . When /? is 
small enough, the system is in a symmetric vacuum, characterised by the formation of FK clusters 
of finite size. If one probes this vacuum with Wilson loops of size much larger than the typical 
dimension of these loops, one finds a perimeter law decay (Wy) °< e^ p ^ [\y\ is the perimeter of 
the loop), the reason being that the only clusters that can be felt by Wy are those near the closed 
path 7 G A. When j3 is larger than a threshold value /3 f which depends on the kind of lattice A, the 
Sq symmetry of the model is spontaneously broken and the corresponding vacuum is characterised 
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Figure 1: A schematic view of the phase diagram of gauge Q-state Potts model. The solid line denotes the 
bulk transition corresponding to the condensation of magnetic monopoles. The dashed line, corresponding 
to the vanishing of the string tension, does not imply any bulk transitions. 



by the formation of an infinite, percolating, FK cluster Goo C G. The spin field a, associated to 
the sites of the lattice A is , in all respects, the disorder parameter of the dual gauge theory and 
the formation of an infinite, percolating, cluster is a direct sign of the condensation of magnetic 
monopoles, i.e. (a,) ^ 0. It is clear that in this case the number of paths of Goo piercing the minimal 
surface encircled by y grows with its area A, therefore one is tempted to argue that large Wilson 
loops obey an area law. Note however that the probe Wy feels only those piercing paths which are 
closed , thus in order to conclude for a confining vacuum one has to assume that also the subgraph 
^ composed by the circuits of G has got an infinite component ^ C "Sf for j8 > j8 f . This has been 



demonstrated through numerical simulations only when Q > 1 [ ]11| , |90. For lesser values this is not 
necessarily true. 

Actually there is a very simple argument showing that, keeping constant the mean number of 
bonds (b) , the size of ^ or, more precisely, the number b<g of bonds belonging to ^ is a decreasing 
function of Q and vanishes in the limit Q — > 0. In fact Q is the cluster fugacity of the system: when 
Q decreases, so does the number of clusters c. The only way to reduce c is to add bridges, i.e. 
bonds that join otherwise disconnected clusters. Now the total number of bonds b is the sum of beg 
bonds belonging to ^ (i.e. to circuits) and of b 0jS bridges, i.e. b = b^ + b 0jS , where G = c € ^ 
therefore a growth of the bridges keeping b constant implies decreasing of beg, q.e.d. 

We note, as a side remark, that there is a general relationship between the bonds of the two 
kinds: 

(b^)^ + (b v )—=N; (3.1) 



this is true for any Q-state Potts model on an arbitrary graph with N links []r2[]. In the special case 
of two-dimensional, infinite square lattice at the transition point (v = y/Q) the self-duality of the 
model requires (b) = ^, thus at criticality we have beg = ^ 2(</Q+i) ' ^" n * s exact result can be used 
as a check of the Monte Carlo algorithm described in the next Section. 

The reduction of closed paths as Q decreases suggests, for Q and /3 — $ small enough, that 
the subset ^ does not longer percolate even in the phase where the symmetry is spontaneously 
broken. For instance, in the limit Q,v — > with the ratio w = v/Q held fixed [ jl3"| ] the surviving 
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configurations are spanning subgraphs not containing any circuits, i.e. ^€ = 0, hence in a 3D lattice 
Wy(G) = 1 , V 7 and G and the Q = dual gauge theory is trivial. 

The above remarks suggest that when Q < 1 is small enough, the standard non-confining vac- 
uum, corresponding to the symmetric phase of the Potts model, is separated from a truly confining 
vacuum, where the subgraph of circuits, holds an infinite component, by an intermediate vacuum 
characterised by the condensation of magnetic monopoles (hence by the formation of an infinite, 
percolating FK cluster) which however is not confining, because the closed paths, corresponding to 
center vortices, form a dilute gas of loops embedded in the infinite cluster rather than a connected 
skein. A sketch of the expected phase diagram of the 3D Potts model in the small Q region is drawn 
in Fig.l. 



4. A Monte Carlo algorithm for Q<1 Potts models 



The non-local cluster algorithm of Swentsen and Wang [14] and its generalisation to non 



integer Q [15 ] is applicable only for Q > 1. In order to study the region < Q < 1 we are interested 



in we have to resort to some local Monte Carlo algorithm Q16| , |17Q . I describe here a variant of the 



method described in Q17[ ] which can be implemented in an efficient way and works only when 
0<<2<l. 

First, divide the interval [0, 1] in three parts a , b — a , 1 — b, where a and b > a > are suitable 
functions of j3 and Q, to be determined later. Then apply the following recursive procedure that 
generates a Markov sequence of spanning subgraphs • • • — * G^ 1 ' — > G^ n+ ^ —* ... of an arbitrary 
lattice A: 

i) Pick a link £ G A and draw a uniformly distributed random number < re < 1 ; 

ii) create or erase a bond on the link £ according to the following rules: (a) put a bond if re < a; 
(b) erase any bond if re > b; (c) in the remaining case (a < re < b) put a bond only if it is a bridge, 
i. e. only if it connects two otherwise disjoint clusters. 

The Markov chain generated in this way forms an ergodic trajectory in the space of con- 
figurations of the Q-state Potts model. Requiring detailed balance with respect the equilibrium 
distribution yields a = 1 — = p and b = The inequality a < b implies Q < 1 1 . 

This Monte Carlo method with its variants was already used to locate the marginal value of 



Q > 2 in three dimensions [ ]17[ ] and in the study of the backbone exponent of critical Q-state Potts 
models in two dimensions [18]. 

Since this kind of algorithms implies a random sequence of disordering moves of type (a) 
and (b), randomly distributed over the lattice, it led to conjecture that they do not suffer of critical 



slowing down Qlq , |17| ]. A subsequent numerical analysis for some integer values of Q in two and 
three dimensions showed that this conjecture is false jl9t]: it reduces critical slowing down, but 
does not completely eliminate it, in the sense that its dynamical critical exponent is smaller than 
Swendsen-Wang, but in general does not vanish, at least for Q > 2 in 2D and Q > 2 in 3D [|TJ]. 

Of course, for a practical use of viable Monte Carlo algorithms it matters not only the intrinsic 
dynamics of the transition rates, but also the efficiency of the numerical implementation. In the 



'if one replaces the rule (c) with the new rule (c'): put a bond only if the number c of clusters is kept constant (this 
was the rule chosen in ref.[M7(] in the case Q > 1) then one finds a = p/Q and b — p. 
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Figure 2: The mass of the lowest physical state in the symmetric vacuum (left) and the string tension in the 
confining vacuum(right). The data refer to Q = 1/10 Potts model in a lattice of size 32 3 . 

present case we succeeded in simulating lattice sizes of the order of those used in current gauge 
simulations. 



5. Numerical results 

A simple observable which can be used to locate the threshold p t = 1 — e~$< where an infinite 
FK cluster forms is the connectivity correlator W(x,y) = (<j>xy), where <p xy = 1 only if x and y 
belong to the same cluster, otherwise is set to zero. Clearly for p < p t one observes an asymptotic 
exponential decay &(x,y) ~ e m x y with increasing separation |jc — y\, where m is the mass of the 
lowest physical state. The correlation length % = 1 jm is of the order of the mean linear size of the 
FK clusters. As a consequence, m vanishes at p t , and it is expected to obey a critical power law 



m ■ 



:a(p t -p) v +b(p t -p) v +... 



(5.1) 



where v = v(Q) is the thermal exponent. Unfortunately it appears that this critical exponent has 
not yet been calculated in 3D Q-state Potts models in the range < Q < 1 , apart the special Q — ► 



j_ 
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We extracted the mass of the 



limit (see the second paper of [|13|]). 

We performed our simulations on a 32 3 cubic lattice at Q 
lowest physical state using the zero momentum projection and fitted the data to (|5J]) as shown in 
Fig.2 (left). As a result the threshold value for the formation of a percolating FK cluster and the 
thermal critical exponent turn out to be p t = 0.0500(18) and v = 2.50(9). Note that this value of 
V is much larger than the corresponding value at <2 = 1 Vq=i = 0.874(2). This agree with the fact 
that in two space dimensions the presumed exact value of v increases as Q decreases. 

On the same lattice at the same value of Q = -L , but at larger values of p, where one can easily 
observe percolation of the sub-clusters made with the circuits of the FK clusters, we measured the 
vacuum expectation value of a set of square Wilson loops, in order to evaluate the string tension. 
They perfectly fitted the expected asymptotic functional form for the confining phase (including 
the log term due to the quantum fluctuations of the underlying confining string). The extracted 
string tension as a function of p is nicely described by a power law (see Fig.2 (right): the fitting 
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curve is a straight line in the log log scale). However the vanishing point of the string tension, 
p„ = 0.063(1), where the deconrining phase starts, does not coincide with the threshold p t . In 
the range p t < p < p the vacuum of this theory is characterised by a non-vanishing magnetic 
monopole condensate which is not confining, being a = there. 

Note that, contrarily to what happens in the gauge duals of Q > 1 Potts models, the critical 
index v a associated to the vanishing of the string tension is totally different from the thermal critical 
index v. It is worth observing that there is no local order parameter that can signal when an infinite 
cluster of circuits forms , being a phenomenon of topological nature that can be detected only by 
large Wilson loops. As a consequence we do not expect that in these gauge models the vanishing 
of the string tension is associated to any sort of bulk transition. 
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